Layer charge instability in unbalanced bilayer systems in the quantum Hall regime 
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Measurements in GaAs hole bilayers with unequal layer densities reveal a pronounced magneto- 
resistance hysteresis at the magnetic field positions where either the majority or minority layer is at 
Landau level filling factor one. At a fixed field in the hysteretic regions, the resistance exhibits an 
unusual time dependence, consisting of random, bidirectional jumps followed by slow relaxations. 
These anomalies are apparently caused by instabilities in the charge distribution of the two layers. 

PACS numbers: 73.50.-h, 71.70.Ej, 73.43.Qt 
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Hysteretic phenomena are widespread in nature. They 
are common magnetic materials, and often indicate a 
non-equilibrium situation associated with a phase transi- 
tion and the presence of domains [l| . Recently, hysteresis 
has also been reported in various two-dimensional (2D) 
carrier systems in semiconductor structures at low tem- 
peratures and high magnetic fields 0, IE El IE IE ■ In 
these cases, magneto- resistance (p xx ) hysteresis appears 
in the quantum Hall (QH) regime when two Landau lev- 
els (LLs) with opposite spin are brought into coincidence. 
While the 2D systems studied have been notably dif- 
ferent, the common thread in these experiments is that 
there is a magnetic transition involving the carrier spin 
1. 

Here we present hysteretic p xx data in 2D bilayer sys- 
tems in the QH regime. The hysteresis in these systems 
has a different origin and is caused by a non-equilibrium 
charge distribution in the two layers. We studied the 
magneto-transport coefficients of GaAs bilayer hole sys- 
tems with unequal layer densities. When the interlayer 
tunneling is sufficiently small, p xx of the bilayer system 
exhibits a pronounced hysteresis at perpendicular mag- 
netic field (B) positions close to where either the majority 
or minority layer is at LL filling factor one. Most remark- 
able is the time dependence of p xx at a fixed field in the 
hysteretic regime, when the two layers are closely spaced. 
As a function of time, p xx exhibits large, random, sudden 
jumps toward higher and lower values, followed by a slow 
decay in the opposite direction. The data may signal an 
instability in the charge distribution of the two layers, 
i.e., an instability associated with the pseudospin (layer), 
rather than spin, degree of freedom. 

We studied nine GaAs bilayer hole samples from six 
different wafers, all grown on GaAs (311)A substrates 
and modulation doped with Si. In all samples, the holes 
are confined to two 15nm-wide GaAs quantum wells 
which are separated by AlAs or AlAs/AlGaAs barriers 
with thickness 7.5 < W < 200nm. The rather thick bar- 
rier combined with the large effective mass of GaAs 2D 
holes [2( reduces considerably the tunneling between the 
two layers [1(1 ]. As grown, the samples have layer densi- 
ties of < 7 x 10 10 cm~ 2 , and low temperature (T) mobil- 
ities of ~ 35 m 2 /Vs. Metallic top and bottom gates were 
added to control the densities in the layers. We stud- 
ied several types of devices, including 2. 5x2. 5mm square 



samples and ones with patterned Hall bars; in these sam- 
ples the ohmic contacts contact both layers. One sam- 
ple was fabricated using a selective depletion scheme 
that allows probing the transport characteristics of in- 
dividual layers. The measurements were performed in a 
dilution refrigerator down to T = 20mK. 

Data of Fig. 1 highlight some of the results of our 
study. In (a) we show a set of traces where p xx was mea- 
sured for a sample with W — llnm as B was ramped up 
or down. For the traces of Fig. 1(a) the total bilayer den- 
sity (ptot) is kept constant at 5.5 x 10 10 cm -2 while charge 
is transferred from one layer to another using back- and 
front-gate biases. We define the charge transfer from one 
layer to another as Ap — (pb — p<r)/2, where pb and px 
are the densities of bottom and top layers, respectively 
fl2|. At a given value of B, we define the filling factor, v, 
of the bilayer system as the ratio between p to t and the LL 
degeneracy, eB/h. We also introduce the filling factors 
for top and bottom layers, vt and vb respectively, as the 
ratio between the layer density and eB/h. 

The data of Fig. 1(a) show that when the bilayer 
system is balanced (Ap=0; top trace) p xx is indepen- 
dent of the direction B is ramped. However, as soon as 
the system is imbalanced (|Ap| > 0) a strong hysteresis 
develops in p xx . For values of Ap < 4.7 x 10 9 cm~ 2 , 
p xx displays hysteretic behavior in two field ranges near 
v = 2, one near vt = 1 and another near vb = 1. When 
Ap > 4.7 x 10 9 cm~ 2 , the hysteresis exists only near 
Vb = 1. For sufficiently large Ap, no hysteresis is ob- 
served. The amplitude of the hysteresis also decreases as 
T is increased (data not shown) and vanishes completely 
above T ~230mK, roughly independent of Ap. 

To probe the contribution of the spin degree of freedom 
to the hysteresis in our bilayer systems, we performed 
measurements in tilted magnetic fields on a sample very 
similar to the one shown in Fig. 1(a). In this experiment 
the direction of the field was kept at an angle 6 with re- 
spect to the normal to the plane of the bilayer system. 
For 8 ranging from 0° to 80° , corresponding to a six- fold 
increase of the total field (and therefore of the Zeeman 
energy) in the hysteretic region, the position of the hys- 
teresis in perpendicular magnetic field did not change at 
all. If the hysteresis were caused by an instability associ- 
ated with the spin degree of freedom, one would expect 
[E that the applied parallel field would change the po- 
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FIG. 1: (a) p xx vs. B traces for a GaAs bilayer hole sample with a barrier width W = llnm. Data are shown for different 
values of charge transfer, Ap, while ptot is kept constant at 5.5 x 10 10 cm -2 . The traces are shifted vertically for clarity. The 
dotted lines indicate the positions of bilayer filling factors v = 2 and v = 1. A strong hysteresis develops in p xx when the 
bilayer system is imbalanced. The right and left tick marks in each trace indicate the estimated positions of filling factor one 
for layers with higher (bottom) and lower (top) densities, respectively. The right panel shows data for a similar sample but 
with independent layer contacts, and p to t = 10.8 x 10 10 cm -2 and Ap = 1.3 x 10 10 cm -2 , (b) p xx vs. B traces measured when 
the ohmic contacts are connected to both layers, (c) and (d) p xx of the top and bottom layers measured individually. In all 
figures the black (red) line represents the trace taken when B is swept up (down). 



sition and magnitude of the hysteresis. These results rule 
out spin as being responsible for the hysteresis. 

To better understand the origin of the observed hys- 
teresis, we fabricated another sample from a different 
wafer, also with W = llnm, using a selective depletion 
scheme |llj | , and aimed to independently probe each layer 
of the bilayer system. The data are shown in the right 
panel of Fig. 1. In (b) we plot p xx for the bilayer system, 
that is, when the ohmic contacts are connected to both 
layers, for both up and down i3-sweeps. In (c) and (d) 
we show p xx traces for the top and bottom layers, mea- 
sured separately, but at the same pair of layer densities 
as in panel (b) traces. Two features of these data are 
noteworthy. First, the traces of (b) exhibit hysteresis in 
two ranges of B. The hysteresis between 1.45 and 1.8T 
matches well the position of z/y = 1 QH state of the top 
layer as seen in (c), while the hysteresis located between 
2.5 and 3.2T overlaps the QH state of the bottom layer 
(see (d)). This observation confirms that the hysteresis 
in p xx of the bilayer system takes place when one of the 
layers is at filling factor one. Second, data of (c) and (d) 
show that each individual layer exhibits hysteresis when 
the other layer is at filling factor one, i.e., p xx of the top 



layer exhibits hysteresis when vg = 1, and vise versa. 

The observation of hysteretic magneto-resistance in 
unbalanced bilayer systems has precedence. Zhu et al. 
[l3j reported hysteresis in 2D electron systems with a 
parallel conducting layer. In their case, the parallel layer 
was a parasitic, low-mobility, doping layer at a distance 
of 75 to 95 nm away from the high- mobility 2D electrons. 
Similar to our data, they observed hysteresis in p xx when 
the layer containing the high-mobility 2D electrons was 
in a QH state. They also presented a simple model to 
explain the observed hysteresis. In their model, the hys- 
teresis comes about because of a non-equilibrium charge 
distribution in the layers. As B is swept, thanks to the 
Landau quantization, the Fermi levels of both layers os- 
cillate. These oscillations lead to temporary imbalances 
between the chemical potentials of the two layers. The 
potential imbalance is particularly abrupt and large when 
the high-mobility 2D electrons enter a QH state as their 
Fermi level jumps by a significant amount, equal to the 
separation between the adjacent LLs. With increasing 
time, of course, the Fermi levels of the two layers have 
to come to equilibrium since the latter are shorted to- 
gether via the ohmic contacts. But this equilibration can 
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FIG. 2: (a) p xx vs. B traces for another sample with 
W = llnm, with p to t and Ap as indicated, (b) p xx vs. 
time measured following an up-sweep of B and stopping B at 
2.195T. p xx displays sudden jumps followed by a slow relax- 
ation. The red lines represent double exponential fits to the 
slow relaxation component of the data, with time constants 
ti ~ 5s and T2 ~ 90s. 



take a long time in the QH effect regime: it has to take 
place via the layers' edges and the ohmic contacts since 
the (bulk) states in the center of the 2D layer that is in 
the QH state are localized and the layer sheet conduc- 
tivity is very small. As a result, p xx , which is recorded 
as B is swept at a finite rate, can show a hysteretic be- 
havior. Consistent with their model, Zhu et al. found 
that when the -B-sweep is interrupted in the hysteretic 
region and p xx is monitored as a function of time, it de- 
cays approximately exponentially toward an equilibrium 
value. Moreover, they found the time constant of the de- 
cay in reasonable agreement with estimates based on the 
parameters of their experiment. 

The hysteresis in our samples resembles what Zhu et 
al. observe and likely has a similar origin. Our data of 
Fig. 1 in fact explicitly show that the hysteresis hap- 
pens when one of the layers is in a QH state, and it is 
the resistivity of the layer which is not in the QH state 
that is hysteretic. The interpretation that the hysteresis 
indicates a charge transfer between the two layers is also 
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FIG. 3: Time evolution of p xx in the hysteretic region for 
samples with different W, ptot, and Ap, as indicated. The 
insets show p xx vs. B for up- and down-sweeps. The arrow in 
each inset indicates the position at which the time evolution 
was recorded after an up-sweep of B. 



supported by our data: the hysteresis observed in p xx of 
either the top or the bottom layer measured separately 
appears much like a horizontal shift of the trace in B, as 
if the layer had slightly different density when ramping 
B up or down [see, e.g., Fig. 1(c)]. The time evolution 
of p xx in our samples, however, is qualitatively differ- 
ent from the observations of Zhu et al., and points to a 
very unusual relaxation. In the remainder of the paper, 
we describe this evolution and speculate on its possible 
origin. 

We have studied a number of samples with varying ptot, 
Ap, and barrier width 7.5 < W < 200nm. Samples with 
W = 7.5nm do not show any hysteresis; this is likely 
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because the interlayer tunneling is sufficiently large so 
that the two layers stay in equilibrium during the B- 
sweep. We observe hysteresis for all samples with 11 < 
W < 200nm, but the time evolution of p xx critically 
depends on W. Examples are shown in Figs. 2 and 3. For 
samples with W = llnm (Fig. 2), the time evolution is 
simply wild! It displays sudden jumps in p xx , followed by 
a slow relaxation after each jump. Note that p xx jumps 
toward both higher and lower values, reminiscent of a 
bistability, although the jumps do not happen between 
fixed values of p xx . In between jumps, p xx follows a slow 
relaxation, in the opposite direction of the jump, that 
can be fitted well by a double exponential [see Fig. 2(b) 
inset]. It is noteworthy that even when measured over 
days (up to 2.5xl0 5 s), we did not observe any tendency 
towards a settling of the jumps. We wish to emphasize 
that, outside the hysteretic region, p xx is independent of 
time to within less than 0.3%. 

We have attempted to quantify the characteristics of 
this evolution by two parameters: average frequency and 
amplitude of the jumps. Our T-dependence measure- 
ments show that, at a fixed B, the average jump am- 
plitude decreases as T increases. Not surprisingly, the 
jumps are no longer visible above the temperature where 
the hysteresis vanishes. On the other hand, the aver- 
age frequency appears to be independent of T. Also, the 
jump frequency and amplitude are independent of the 
magnitude of the sample current, as long as the current 
is kept sufficiently small (< 10 nA) [lj|. Interestingly, 
the jumps and decays appear to continue even when the 
current is completely turned off |l5j . 

The sample with W = 30nm shows a behavior qual- 
itatively similar to the one with W — llnm, although 
both the frequency and size of p xx jumps are smaller 
[Fig. 3(a)]. Data for the W = 50nm sample, however, 
are qualitatively different [(Fig. 3(b)]: there are by far 
fewer p xx jumps (typically one jump every few 10 3 s), 
and p xx appears to decay with time. Finally, for the 
W = 200nm sample [Fig. 3(c)], we typically observe a 



simple decay with time |l6j . Interestingly, for the sam- 
ples of Figs. 3(b) and (c), p xx continues to decay with 
an ever increasing time constant. 

The time evolutions we observe for the sample with 
W = 200nm barrier [Figs. 3(c)] is qualitatively similar 
to the observation of Zhu et al. 13]. The time depen- 
dences for the samples with barrier widths W = 11 and 
30nm, however, are very unusual and cannot be under- 
stood in a simple model where the bilayer system slowly 
and steadily relaxes to an equilibrium state. In these 
samples, p xx displays sudden jumps that do not have any 
tendency to settle, at least over a time scale of days. We 
do not know the origin of these time evolutions. The sud- 
den jumps in p xx in our bilayers with small W bear some 
resemblance to the so-called Barkhausen jumps, which 
are observed in magnetic materials pj. The Barkhausen 
jumps occur when the magnetic system finds a lower en- 
ergy state available and one or several domains change 
orientation. A tantalizing speculation is that in the bi- 
layer systems with close layer separation, the interlayer 
interaction acts as an opposing force to the charge trans- 
fer caused by the Fermi level difference ■ I n this 
scenario, the two opposing mechanisms may mediate the 
creation of a complicated layer charge density pattern or 
pseudospin domains. 

Another possibility is that the observed jumps are not 
intrinsic to the sample, but rather triggered by external 
sources (e.g. electromagnetic noise). If so, it is a puzzle 
why the jumps are much more frequent in the bilayer 
samples with smaller layer separation and are seen only 
in the hysteretic region. 

In summary we report an unusual time dependence 
associated with hysteretic magneto-resistance in GaAs 
bilayer holes with close layer separation. The resisitivity 
exhibits sudden jumps with time, possibly caused by a 
layer charge instability. 
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